Glueball-quarkonia content and decay of scalar- isoscalar mesons 



< 

(N 



M. Strohmeier-Presicek^ , T. Gutsche^, R. Vinh Mau^ and Amand Faessler^ 

^ Institut fiir Theoretische Physik, Universitdt Tubingen, Auf der Morgenstelle I4, 

D-72076 Tubingen, Germany 

d> , „ 

Q> ' Laboratoire de Physique Theorique des Particules Elementaires, 

Universite P. et M. Curie, 4 Place Jussieu, 75252 Paris Cedex 05, France 

Abstract 

We investigate the hadronic two-body decay modes of the scalar-isoscalar 
/o(1370), /o(1500) and /o(1710) states as resulting from the mixture of the 
^ . lowest lying scalar glueball with the isoscalar states of the ground state "^Pq 

^ : 

■ QQ nonet. In the decay analysis we take into account the direct coupling 

. of the quarkonia and glueball components of the /o states to the quarkonia 

components of the two-meson final state with the decay dynamics inspired 

Oh! 

D . by the strong coupling limit of QCD. We calculate partial decay widths for 

the /o states in the proposed three-state mixing schemes and discuss their 
compatibility with the observed decay features. Finally, we determine the 
glueball-quarkonia content of the /o states from a detailed fit to experimental 
decay data of /o(1500) and give predictions for the partial decay widths of 
/o(1370) and /o(1710), providing thus a sensitive test of the proposed mixing 
scheme. 
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I. INTRODUCTION 



The identification of mesons beyond tlie conventional valence quark- ant iquark {QQ) pair 
configuration is one of the most challenging tasks in low-energy meson spectroscopy. 
For instance, as a consequence of gluon self-couplings in QCD the usual QQ spectrum is 
expected to be supplemented by mesonic states with pure gluonic composition, that is glue- 
balls 0]. 

Lattice QCD, in the quenched approximation, predicts the lightest glueball state (Go) to 
be a scalar {J^^ = 0++) lying in the mass range of 1.4 - 1.8 GeV [0,^. Since the observed 
isovector and isodoublet states ao(1450) and Ko(1430) set the natural mass scale of the 
ground state scalar meson nonet, the isoscalar QQ partners nn = -^{uu + dd) and ss are 
expected to lie in a mass region close to the predicted ground state glueball. The vicinity of 
masses for the pure Gq, nn and ss mesonic states suggests that the intrusion of the scalar 
glueball into the scalar QQ meson spectrum is accompanied by significant mixing of Gq with 
the isoscalar members of the ^Pq nonet 

Several glueball-quarkonia mixing schemes have been proposed to account for this 

scenario and to reveal the Gq/QQ nature of the scalar-isoscalar states /o(1370), /o(1500), 
established by Crystal Barrel at LEAR 0] in proton-antiproton annihilation reactions, and 
the /o(1710) 0, which is the scalar component of the /,=o,2(1710) |T^. In a simplified 
picture considered in Refs. the scalar glueball ground state Go mixes with the lowest 
lying ^Pq isoscalar quarkonium states. Further admixture of the first radially excited ^Pq 
nonet is introduced in Refs. [^0. However, the restriction to the glueball and near lying 
^Pq quarkonia ground states is a good approximation provided that the excited glueball and 
quarkonia states are high enough in mass. 

Quantitative predictions in the three-state mixing schemes of Refs. 0,^] for the glueball- 
quarkonia content of /o(1370), /o(1500) and /o(1710) differ substantially due to orthogonal 
theoretical assumptions concerning the level ordering of the bare states before mixing. While 
Ref. m uses for the input masses of the bare states M{nn) < M{Gq) < M(ss), in Ref. 
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the pure glueball lies above the scalar QQ state. The difference of level ordering in the 
input bare masses leads to a substantially different Gq/QQ content of the physical /q states, 
especially for /o(1500) and /o(1710). While in Ref. the glueball is distributed with nearly 
equal strength over the three /o states and /o(1710) has a large ss component, the analysis 
of Ref. leads to a dominant Gq structure in the /o(1710), while /o(1500) is a nearly pure 
ss state. The results for the /o(1370) are comparable in both approaches, with a large nn 
component residing in this state. The observed total width of the /o(1370) ^ and its strong 
coupling to the vrvr decay channel is in good agreement with naive quark model expectations 
|[TT| , P^ for a nn ^Pq state and seems to confirm its dominant nn structure. Therefore, the 
main difference between the two mixing schemes is the amount of glueball and ss strength 
residing in the /o(1500) and /o(1710) states, respectively. A further test of the proposed 
mixing schemes, where the glueball ground state intrudes in the scalar quarkonium sector, 
is the analysis of the decay modes of the observed /o states. The decay of /o(1500) into two 
pseudoscalar mesons was studied in leading order of strong coupling QCD Neglecting 
possible gluonic components of the isoscalar mesons in the final state, the vrvr, KK, rjr] and 
ry?7'decay modes are entirely driven by the quarkonia components of the /o(1500) state. In 
this approximation, results for the partial decay widths deduced from SU (3) flavor sym- 
metric couplings are quite similar for both mixing schemes and cannot be discriminated , 
e.g. the observed weak KK decay mode is obtained in both schemes by the destructive 
interference between the nn and ss components. 

In order to reveal the difference of the proposed three state mixing schemes, one should take 
into account the direct coupling of the Gq component to the final state mesons. This mech- 
anism occurs in the leading of strong coupling QCD by the transition Gq — > GqGq, which 
can produce final states with isoscalar mesons {rji], rji]' and aa, where a is the broad low lying 
TTTT S-wave resonance) with possibly a sizable gluonic component. Independent of a possible 
small glueball component in the final state isoscalars is the decay process Go {QQ){QQ)- 
This next-to-leading order transition can contribute to all two-body decay modes of the 



mixed /o states. An investigation |jT3| of the two-body decay modes aa, pp and 7r*(1300)7r 



of the /o(1500), all leading to the An decay channel, indicates the possible presence of a 
sizable direct coupling of the gluonic component of the /o(1500) to the An decay channels. 
Neglecting a possible Go component in a, the leading order decay mechanism predicts 
the hierarchy of branching ratios BR with BR{pp)>BR{-K'K)>BR{aa) > fii?(7r*(1300)7r), 
independent of any particular three-state mixing scheme. The results are in strong conflict 
with current Crystal Barrel data, where the pp decay mode of the /o(1500) is weak or even 
absent JlSl and a stronger coupling to the An modes than to the 2tt decay channel |^ is 



required. 

In the present work, going beyond the lowest order decay mechanism, we couple the QQ 
components of the mixed /o states to the quarkonia components of the final state mesons. In 
the decay analysis we also take into account the transition Go (QQ) (QQ) occurring as a 
next-to-leading order decay mechanism deduced from strong coupling QCD p. Inclusion of 
the direct coupling of the Go component of the mixed fo states to the quarkonia content of 
the final state mesons will lead to contradictory physical consequences in the proposed three- 
state mixing schemes of Refs. [^,0, in particular for the physical state /o(1710). Because 
of the shortcomings of these proposed mixing models, we undertake to extract a three-state 
mixing scheme based on a fit to the detailed experimental data of /o(1500) two-body decays 
and give the predictions for partial decay widths of the partner states /o(1370) and /o(1710). 
This paper is organised as follows. In Sec. |I| we introduce the three-state mixing scheme 
of the lowest lying scalar glueball with the quarkonia states and the decay dynamics as 



suggested by the strong coupling limit of QCD. In Sec. |T| we develop the formalism for 
the decay of the /o states by its QQ components in the framework of the pair creation or 
^Pq model. In Sec. |I^ we construct a digluonium wave function based on cavity QCD for 
the effective description of the scalar glueball state. Furthermore, we indicate the transition 
mechanism for the decay of the Go component into (QQ) (QQ)- The results are presented 
in Sec. |V|; we first show how distinct decay patterns result from the unmixed QQ and Go 
components of the /o states. We then give partial widths for the two-body decay modes of 
the fo states in the mixing schemes of Refs. and confront the results with experiment. 



Finally, we present a three-state mixing scheme deduced from a detailed fit to the experi- 
mental data for /o(1500) decays and give predictions for the partial decay widths of /o(1370) 
and /o(1710), providing thus a sensitive test of the proposed mixing scheme. A summary 
and conclusions are given in Sec. |V|. 



5 



II. THREE-STATE MIXING AND DECAY OF SCALAR-ISOSCALAR STATES IN 

STRONG COUPLING QCD 



In the lattice formulation of the strong coupling limit {g oo) of QCD [^,n|, glueballs 
and QQ mesons are noninteracting eigenstates of the QCD Hamiltonian. Departure from 
this limit causes mixing and decay of the strong coupling eigenstates 

To lowest order, the mixing of the scalar glueball Go and the ^Pq quarkonia states nn = 
1/ J (2) (uu + dd) and ss can be introduced by the interaction Hamiltonian Hi 



Hi 



z niss 
\/2z m„i 



(1) 



where z =< Go\Hi\ss >=< Go\Hi\nn > /v^ represents the flavor independent mixing 
strength between the glueball and quarkonium states. The masses of the bare states i = 
GQ,ss,nn before mixing are denoted by rrii. Deviation from the SU{3) flavor symmetry is 
known to be small |T^ and quarkonia mixing is assumed to be a higher order perturbation 
in the strong coupling eigenstates |Go >, \ss > and \nn > All parameters of Hi can be 
taken to be real and positive. 

Hi possesses three eigenstates |/o(M) > with physical masses M , which are are given by 
the linear combinations: 



|/o(M) >= ann\nn > +ass\ss > +aco\Go >, 



(2) 



with a^^ + + a^^ = 1. Choosing the bare mass to be equal to that of the ob- 
served isovector state ao(1450) |10[, the remaining parameters in Hi can be adjusted to 
give the masses of the observed resonances /o(1370), /o(1500) and /o(1710), which flxes 
the quarkonia-glueball composition of the physical states 0. In Ref. [^] it is assumed that 
the level spacings of the bare states fulfll the relation {tjigo — uinn)/'^ = ^ss — "^Gq = 
which again determines the state vector of the physical fo{M) states. In this case /o(1500) 
is degenerate in mass with the Go state and the physical masses M of the mixing partners 



/o(1370) and /o(1710) depend on the mixing strength z. 

However, there are sizable uncertainties in these procedures, particularly since the mass of 



the /o(1370) is poorly determined, ranging from 1200 - 1500 MeV |T0[ and the exact level 
spacing of the bare states in unknown. For these reasons we will also present an alternative 
strategy by assuming that /o(1500) is one of the eigenstates of Hi and by extracting the 
mixing coefficients {i = nh, ss, Go) directly from the experimental two-body decay data. 
This will allow us to give the bare masses rrii entering in Hj of Eq. (1), in terms of the 
mixing strength z, which in turn is restricted by the physical masses of the partners of the 
/o(1500). The resulting glueball-quarkonia composition of these partner states , however, is 
completely independent of z, solely determined by the mixing coefficients of the /o(1500). 
In the minimal constituent picture Gq consists of two gluons bound in a color singlet con- 
figuration with J^'-' = 0^~^. For finite g the mixed |/o(M) > decays in lowest order by its 
Go and QQ components. The relevant diagrams are shown in Figs. ^ and ^ The transition 
Go Go Go of Fig. 1 couples the Gq component of the physical state |/o(M) > to the Gq 
component of the final state isoscalar mesons. The constituent gluons of the initial state 
decay into two pairs of low energy gluons which involves two three-gluon vertices in the 
decay mechanism. We will neglect this transition in our decay analysis by assuming that 
the considered final state isoscalar mesons do not contain a Gq component. In principle, 
the transition of Fig. 1 can contribute to the decay modes /o(M) — > rjrj^rjri' and crcr, de- 
pending on the unknown overlap of the rj, rj and a with the corresponding glueball states. 
Current evidence for the relevance of such a process involving the pseudoscalars ?7 and rj is 
not compelling, as for example discussed in Ref. [Q. For the scalar a meson a large coupling 
to the glueball state is not excluded @,|10|- Therefore, we should keep in mind that the 



Go Gq Gq transition could induce some changes for the r]r], rjrj' and aa decay modes of 
the /o states. 

The leading order decay mechanism QQ {QQ) (QQ), illustrated in Fig. H, leads to the 
decay of |/o(Af) > by its quarkonia components QQ and is familiar from the OZI-allowed 



meson decay |Tl|J13||r6[| . This process will be described in Sec. |ITl| in the framework of the 
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Pq pair creation model [17 . 



Here we also take into account the next to leading order decay mechanism Go (QQ) (QQ) 
as indicated in Fig. ^ In this mechanism the Go component of the physical fo{M) states 
couples to the final state QQ mesons by conversion of the constituent gluons into quark- 
antiquark pairs. Sec. fV| is devoted to the description of this process. 

Neglecting a possible Gq component in the final state isoscalar mesons, the partial decay 
width Tfg^BC for the decay of \fo{M) > into a two-meson state B C can be written as 

EaEr 



where M^q^^^ and Mq^^^bc denote the transition amplitudes including phase space of the 
respective decay mechanism. The sum extends over the relative orbital angular momentum 



Ibc between the mesons B and C. With the decay momentum \K\ = K, Ei = y Mf + K"^ 
is the energy of the final state mesons i = B,C with mass Mj. The strengths of the quark- 
gluon coupling and the pair creation amplitude in the ^Pq model are factorised out and 
denoted by 7 and A, respectively. The relative phase between the two transition amplitudes 
is fixed by choosing M^q^^^ and Mq^^^J^^j to be real and the ratio 7^/ A to be complex with 

^ = Ke't (4) 
A 

In this case Tf^^BC is given by 

T/o^BC = |m^^^^^ + M^^^^^ + 2 K cos<P M^'l^^^^M^j--),^! . (5) 

with the summation over Ibc included in the definition of Mqqi^g^^^^bc- Fmally, in order 
to take into account properly the available phase space for the decay /o — > BC we average 
over the mass spectrum /(m) of broad mesons: 

Tf^^BC = J drnf^drriBdmc T f.^^Bci'mf^, tub, mc)f{mfg)f (1713) f{mc) (6) 

(r./2)^ 



/(m) oc 



with a threshold cutoff as in Ref. |T^. The individual masses Mi and widths Fj of the 
resonances are taken from the Particle Data Group. 

The mass distribution of the scalar-isoscalar a meson is parameterised with Mo- = 760 MeV 
and = 640 MeV [|1|]. For the 7r*(1300) we use the resonance width of T^. (1300) = 200 MeV 
. The relationship between /o(1370) and the broad structure around 1100 MeV H 



called /o(400 — 1200) [|l^ is not entirely clear. If the two states are manifestations of a 
single object, then rjQ(i37o) ~ 700 MeV For two independent states one has rjQ(i37o) = 
351 ± 41 Mer,My„(i37o) = 1360 ± 23 MeV § and 1/0(1370) = 230 ± 15 MeV^, My„(i37o) = 



1300 ± 15 MeV 28 
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III. DECAY OF THE QUARKONIA COMPONENTS 



The decay of the quarkonia components of the |/o(M) > states into the final state mesons 
B and C is calculated in the non-relativistic ^Pq or pair creation model |l7]. This model has 



a solid foundation in strong coupling QCD |rT|PT| , P^ . It describes the OZI-allowed decay 
QQ — ^ BC by the creation of a QQ pair with quantum numbers {I'^{J^^) = 0^(0++)) out 
of the hadronic vacuum, as indicated in Fig. ^ In the evaluation one also has to consider a 
second diagram where the quark 3 (4) goes into the meson C (B). The constituent quarks of 
the initial state are spectators in the transition. The strength of this transition is governed 
by the dimensionless constant A, which is related to the pair creation probability. From 
detailed fits to tensor meson decay p|,p3| it is concluded that the pair creation mechanism 



is flavor independent. 

The nonperturbative QQ ^Pq vertex can be defined as 

^(43) 



where ^4(3) are the momenta of the quark (antiquark) 4 (3) and 3^i(p) = \p\ Yi{p). The 
identity operators If and Ic project onto singlet states in fiavor (F) and color (C) space of 
the created QQ pair (43). 

Using the harmonic oscillator ansatz '^n=o,i=i{'Pi,P2), the wave function of the QQ compo- 
nent of the mixed /o can be expressed in its centre-of-momentum frame as: 



(43) 1(43) 



^/o(Pl,P2) = +P2) [X5=l(12) (g) ^„=o,«=i(pi,P2)]j=o (12) Xc 
"9 f?5 

_3v/?_ 



6{pi + P2) exp [-^RIq{pi - P2f] (8) 



where Rqq is the size parameter for the QQ component. The intrinsic spin X5=i(12) of the 

QQ pair (12) is coupled with relative angular momentum / = 1 to the J^*" = 0^^ scalar 

(12) 

ground state (n = 0) and Xc the color singlet wave function. The SU{3) fiavor part 
Xp{l2) is given by 
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x|?(12) = \n1n2 > + ass \s1S2 >, (9) 

with the quarkonia mixing coefficients aj [i = nn, ss) of the three-state mixing scheme 
introduced in Sec. ||. For the final state mesons i {i = B, C) with momentum Pj, total spin 
Jj, relative orbital angular momentum U and intrinsic spin Si we use the wave functions 

where pi(m) are internal quark momenta for the quarks l{m). Explicit forms for the harmonic 
oscillator wave functions "^mk used can be found in Appendix A. 

The transition amplitude M^^'^„„, as defined by Eq. (3), separates into a fiavor part 



M^Q^sc = <xl(13)x^(42)|lfV^?(12)> (H) 



and the spin-spatial-color part Mqq^2,^bc 



J ^ / n dp. ^(P^ + P3) S{p, + n) (12) 
5{pi +P3-K) 5{p2 +PA + K) [x^Sb ® ^n^iAPi^n)]^ 



(43)t 



,(12) 



00 

with the decay momentum K and the color matrix element < Xc'^^'Xc'^^ I IcIXc^'' ^ ~ \J'^I'^- 
The amplitude M^q^^^ is evaluated analytically. Since there exists exhaustive hterature 
iTl| , [l6|Jl7| on the ^Pq model we just cite the results in Appendix B. Values for the various 
fiavor matrix elements are given in Table |. The results of Table H were obtained by assuming 
X'f — {nfi >. Thereby we have used the following definitions for the fiavor content of the r) 
and T]' mesons: 

ipp = aps \nn > — [3ps \ss > (13) 
= Pps \ns > + aps \ss > . 

The quantities aps and Pps are simply related to the pseudoscalar mixing angle 6'p as defined 
by the Particle Data Group |jTO|. For numerical results we use the value of = — 17.3" 
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g or equivalently aps = 0.794, /3ps = 0.608. 



The partial decay widths for the quarkonia decay QQ BC are defined as in Eq. (5) 
and depend on the pair creation amphtude A and the quarkonia size parameters Rqq and 
Rb = Rc for the decay channels considered here. 

Exhaustive experimental data exist for the hadronic decays of the ground state tensor meson 
nonet and the strategy is to extract the meson size parameters combined with the ^Pq 
strength from a full scale decay analysis. A recent fit to tensor meson decays yields 
= 2.11 and Rb = 3.65 GeV~^ for nonstrange S-wave mesons and strange tensor mesons, 
suggesting Ra^ = Rqq = 3.65 GeV^^. S-wave mesons containing strange quarks have an 
effectively reduced size with Rb = 2.48 GeV~^ for t] and rj' and Rk = 2.82 GeV~^ for K. 
No fiavor symmetry breaking at the ^Pq vertex was obtained, hence A is independent of the 
created QQ fiavor. Due to orthogonality of the vr and 7r*(1300) wave functions the respective 
size parameters fulfil the relation R^, = R.„*. The R^ dependence of the quarkonia decay 



amplitude was shown to be small [|T^. We therefore choose R^j = i^jr for our numerical 
results presented in the forthcoming section. 
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IV. DECAY OF THE GLUEBALL COMPONENT 



In this section we consider the two-body decay of the ground state scalar glueball (Go) 
in the next-to-leading order decay mechanism of strong coupling QCD For the decay 
analysis we adopt a non-relativistic picture where Go is a bound state of two massive con- 
stituent gluons g. The decay is proceeding via the conversion of the constituent gluons into 
QQ pairs as illustrated in Fig. ^. In the evaluation a second diagram with different quark re- 
arrangement in the final state analogous to quarkonia decay has also to be considered. The 
elementary quark-gluon interaction vertex is given to lowest order in the non-relativistic 
limit by 

Vi^^Q.Q.) = ^ S{i - m - Pk) {a^ik) ■ 61) if ^ {\ E (14) 

^ a=l 

with the internal momenta qi for gluon i=l,2 and pn^k) for the quarks with label 1=1,3 (k=2,4). 
The identity operator if ^ projects onto a flavor singlet state of the created QQ pair (Ik). 
The last term in Eq. (15) is the color part of the interaction vertex with the Gell-Mann 
matrices A'^ acting in color space of (Ik). The color octet wave function of the gluon i with 
polarisation vector is denoted by A". The strength of the gi QiQk transition is given 
by 7. We assume that flavor symmetry is manifested in the quark-gluon coupling [^], hence 
< nn\V\g >=< ss\V\g >= 7. In order to evaluate the decay amplitude Mg^^bc we will 
construct a non-relativistic digluonium wave function from cavity QCD. 



A. Digluonium wave function 

In constructing the glueball wave function we use the Bag Model as a guidance. The 
ground state scalar glueball is realized in the Bag Model by the coupling of two transverse 
electric (TE) modes of the constituent gluon {J^*" = 1^~) to scalar quantum numbers 
J^'" = 0^^ 1^^. Using the non-relativistic TE gluon mode, proposed in Ref. ||26|, the scalar 
glueball wave function is in simple correspondence with the Bag Model solution and can be 
written as 
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^Go(n,r2) = [T^i(l+-) ® T^2(l+-)]j=o (15) 
= { [^n,=oh=iin) ® e^]i ® [^„2=o/2=i(r2) ® e^'] ^ }^^^ ^^^^ 

where \l/„=o«=i(^ is the single particle harmonic oscillator wave functions as defined in 
Appendix A. The spherical components of the gluon polarisation vector e are given by e^^'^'\ 



(12) 

The color singlet wave function for the two gluon state is 



where we have chosen the "cartesian" basis A'' with 



V2 ij^^ 



(16) 



(17) 



for the octet representation of color ^[/(S) with Cj and Cj being the basis states for the 
fundamental representations 3 and 3. 
^Go(^i!'^2) can be recoupled as 

1 1 1 



2(7+1 



(18) 



1 1 9 

[^„,=o/i=i(ri) ® ^„2=oz2=i(r2)]3 ® 



^ 



(12) 

C ■ 



^9 J J=0 

Introducing relative and centre-of-mass coordinates r = r i — r2 and R = l/2(ri + r2) 



respectively, we perform the expansion 



[^'ni=o/i=i(n) ® ^„2=o/2=i(r2)] = ^ < n/, NL,g\ni = Oh = 1, 7^2 = O/2 = 1, 5- > 



nlNL 



(19) 



where energy conservation restricts the sum to 2n + / + 2A^ + L = 2. The relevant values for 
the transformation brackets < nl, NL, g\nili,n2l2, 9 > can be found in Appendix C. 
When projecting out the spurious centre-of-mass excitations, the normalised glueball wave 
function can be written in its centre-of-momentum frame as 



^Go('?l,'?2) 



^ 



2 Jo 



(20) 



{ ^n=iz=o(gi, 92) ® [et ® ei]^ }^ j 5{qi + ^2) ^ 



,(12) 

C 
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with the internal gluon momenta (i = 1, 2) and ^2) given in Appendix A. 



B. Digluonium decay amplitudes 

With the digluonium wave function of Eq. (21) were are now in the position to evaluate 
the decay amplitudes for the process Go ^ BC of Fig. |. The decay amplitude M^^Zbc: 
as defined by Eq. (3), separates in analogy to quarkonia decay into a flavor part Mq^^^^j 
with 



<xf(13)x?(42)|ir)ir)|0> 



(21) 



and a spin-spatial-color part 



M, 



{Ibc)SSC 



Go-^BC 



V27tK J f ^ BC I dQidQ2 H ^Pi S{qi - pi - P2) 

V ^^ho J i=i,.A 

S{q2 -Pa-Ps) *c(P4,P2) ^Upi^Ps) (0^(12) ■ e!) (0^(43) ■ £2) ^Go{Qu?2)- 



(22) 



The flavor matrix elements in the SU (3) flavor symmetry limit are given in Table |T[ The 
color transition amplitude Mq^^^q is given by 



m; 



c 



Go^BC — ^<Xc Xc 1X^-^(12)^1 X! -^Ol3)^2l^C > 



d=l 

24v^ ti 3 



The analytic evaluation of Eq. (22) contains the spin matrix element 



(23) 



< X5s(13) X5c(42)| ((T(i2) • el) ((T(43) ■ ei)| 



5=0,2 



> 



6a/(2Sb + 1)(2^c + 1) 



\5c+l 



< SbMbScMc\Sh + fi' > 



1/2 1/2 Sb 
1/2 1/2 5c 
1 1 S 



(24) 



and the spatial overlap 
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/ dqidq2 H ^^Pi - Pi - P2) S(q2 - P4 - Ps) S(qi + 52) (25) 

1=1, ..4 

■S(pi+P3 - K) 5{P4+P2 + K) ^igig{p4,P2) *ic/c(^l'-^3) *n=0,U=2,o(gi,g2) 

= Si^;0 Smffi ^ J dqi dp "^n^i^qi - 1/2{K - p), -qi - 1/2{K + p)) 

where If is the relative orbital angular momentum of the final state with projection m/. A 
complete listing of the analytic results for the amplitudes Mq^^^q is given in Appendix D. 
For the decay Gq BC {If) we obtain the dynamical selection rule that the final state 
partial wave is restricted io If — Hence the decay fo{M) — > oiTT will proceed entirely by 
the quarkonia components of fo{M). Moreover, the D-wave decay amplitude {If — 2) for 
the two vector meson final state is forbidden in contrast to quarkonia decay. 
With the QQ meson size parameters fixed, Mq^^^q depends on the glueball size parameter 
Rqq . In the constituent picture we expect the mass Mq^ of the digluonium to be ~ 2 nig 
with a lattice motivated constituent gluon mass of about rUg ~ 750 — 850 MeV. Assuming 
equal confinement strength for quarks and gluons on the constituent level we get Rgq — 
3.65 /V2^ GeV~^. With this choice the constituent mass for the hghtest quark flavors lies 
in a reasonable mass range of 300 — 340 MeV . 
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V. TWO-BODY DECAYS OF SCALAR-ISOSCALAR STATES AND 

THREE- STATE MIXING 



In the following we discuss the application of the decay models developed in the previous 
sections to the phenomenology of the /o decays. Given the detailed experimental decay data 
on the /o(1500), our analysis will focus on this state with decay predictions given for the 
partner states /o(1370) and /o(1710). After a brief overview of the experimental results for 
/o(1500), we will illustrate the distinct decay patterns resulting from the unmixed QQ and 
Go components. In a first analysis we apply the decay models to the three-state mixing 
schemes constructed in Refs. 0,0. Then we directly extract the mixing coefficients from a 
fit to the experimental two-body decay data. 

A. Observed decay properties of /o(1500) 

The scalar-isoscalar resonance /o(1500) is now clearly established in proton-antiproton 
{pp) annihilation reactions by the Crystal Barrel collaboration at CERN |^ . Mass and width 
of /o(1500) are determined as [|]: 

M/„(i5oo) = 1505 ± 9 MeV T^^dsoo) = HI ± 12 MeV. (26) 

For the two-body decay modes /o(1500) BC Crystal Barrel obtains the branching ratios 
BR{BC): 

BR{nn) = 29 ± 7.5% BR{riri) = 4.6 ± 1.3% BR{r]T]') = 1.2 ± 0.3% 
BR{KK) = 3.5 ± 0.3% BR{A7t) = 61.7 ± 9.6%. (27) 

Here, the An final state includes aa and 7r*(1300)7r as intermediate decay channels. The 
allowed pp decay mode of the /o(1500) is found to be weak ||T3[. Taking the central value for 
the width we deduce from the measured branching ratios the partial decay widths r{BC): 

T{27i) = 32.2 ±8A MeV T{r]r]) = 5.2 ± 1.5 MeV T{r]r]') = 1.4 ± 0.4 MeV 
T{KK) = 3.9 ± 0.4 MeV T{An) = 68.5 ± 10.7 MeV. (28) 

17 



Alternative multi-channel analyses based on data sets from previous experiments have 
also been performed. The results for the KK and rjrj modes are within the experimental 
errors of the Crystal Barrel data. A larger total width {Ttot = 132 ± 15 MeV) and tttt decay 
width (Ftttt = 60 ± 12 MeV) is obtained suggesting a smaller value for r(47r) than indicated 
in Eq. (28). 



B. Two-body decays of the QQ and Go components 

We first discuss the decay patterns resulting from the individual decay mechanisms of 
Figs. and coupling to the QQ and the Gq components of the /o states, respectively. 
In lowest order, neglecting a possible Gq GqGq contribution, the decay /o — * BG will 
proceed via its quarkonia components QQ. In general we consider a mixed quarkonium state 

\QQ > = cosa \nn > — sina \ss > . (29) 

The dependence of the branching ratios BR for the two-body decays of /o(1500) on the 
mixing angle a is given in Fig. ^. The results of Fig. | correspond to the original discussion 
in Ref. for the two pseudoscalar meson decay modes of a ^Pq isoscalar QQ state. But 
now predictions for all possible two-body decay modes are given. 

With the three-state mixing schemes of Ref. {ann = 0.43, a^j = —0.61, corresponding to 
a = 54.8°) and Ref. {ann = 0.40, a^s = —0.90, that is a = 66°) we obtain following ratios 
for the pseudoscalar decay modes 

BR{mT) : BR{KK) : BR{r]r]) : BR{7]7]') 

= 1 : 0.21 : 0.005 : 0.93 Ref. [4] (30) 
= 1 : 0.97 : 0.12 : 1.80 Ref. § (31) 

to be compared with the Crystal Barrel results of 

BRiTTTv) : BR{KK) : BR{'r]r]) : BR{r]r]') (32) 
= 1 : 0.12 ± 0.03 : 0.16 ± 0.06 : 0.04 ± 0.02. 
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Best agreement with data is achieved for /o(1500) ~ \nn >, for which we obtain 

BR{'K^) : BR{KK) : BR{t]T]) : BR{r]ri') = 1 : 0.19 : 0.16 : 0.10. (33) 

Hence, experimental data for the ratios of two pseudoscalar decay modes of the /o(1500) 
are consistent with a dominant nn interpretation. However, the predicted total width for 
a pure ^Pq nn state is in conflict with experimental observation. In Fig. ^ we indicate the 
calculated total width for a /o(1500) state with the QQ configuration of Eq. (^91) . A pure nn 
state results in Ttot ~ 500 MeV, which is considerably higher than the observed width. A 
sizable reduction of the QQ amplitude in the /o(1500), as affected by the admixture of a Go 
component, leads to a reduction of the total width. For example, the mixing schemes of Refs. 
@| and lead to a total width of F = 111 MeV and F = 127 MeV, respectively. Working 
in the lowest order of the decay mechanism, i.e. /o(1500) decays by its QQ components, the 
observed two pseudoscalar decay modes are consistent with a dominant nn content, slightly 
favoring the mixing scheme of Ref. 0]; the possible presence of an admixed Gq component 
is then only reflected in the corresponding reduction of the total width. 
When going beyond the two pseudoscalar meson decay modes in the leading order decay 
scheme, strong deviations from the observed decay pattern of the /o(1500) arise [|T3|. For 



the two-body decay modes leading to the An decay channel, including also charged pion 
combinations, we obtain 

BRi^TiTi) : BR{pp) : BR{7i*n) : BR{aa) : BR{ai7i) = 1 : 0.97 : 0.14 : 0.31 : 1.8. (34) 

Predictions are independent of the mixing angle a, since the final states of Eq. (^) only 
couple to the nn configuration. Relative to the vrvr channel the aivr mode is the strongest 
decay channel. The predicted large pp branching ratio with BR{pp) ^ BR^-k-k) is in conflict 
with the preliminary analysis |]T5| , where the pp channel is found to be strongly suppressed. 
Furthermore, for the combined aa and 7r*7r decay channels with BR^-k* an) =0.1 [^] we 
deduce the ratio BR{4:Tt)/ BR{2tt) = 0.32, which is much lower than the experimental result 
of B R{4:Tt) / B R(2tt) = 2.1 ± 0.6 [§. The observed suppression of the pp and enhancement 
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of the An'^ decay modes are both in strong conflict with a naive QQ interpretation of the 
/o(1500) and give further indication for a sizable direct coupling of the gluonic component 
to the decay channel. 

Next we consider the decay pattern evolving from the direct decay of the Go component as 
indicated in Fig. 0. Branching ratios for the scalar glueball (Go) decay in dependence on 
its mass are shown in Fig. |^. A main signature is the strong suppression of the two vector 
meson decay channels, particularly the pp, over the entire mass range, due to the dynamically 
forbidden D-wave amplitude. The 11*71 and aa decay modes exceed the contribution of the 
TTTT decay channel for a Go mass of about 1500 MeV. These results are in clear contrast to 
those obtained for the QQ decay mechanism. Both, the suppression of the pp decay mode 
and the enhancement of the 47r° decay channels resulting from the decay of the gluonic 
component Go are features observed in the /o(1500) decay pattern as discussed previously. 
For the two pseudoscalar meson decays the results differ from the naive expectation 

BR{7i7i) : BR{KK) : 5i?(r/77) : BR{7]ri') = 1:^:^:0 (35) 

deduced from flavor symmetric couplings. For a glueball mass of M^-y = 1500 (1700) MeV 
we obtain 

BR{'KTi) : BR{KK) : BR{r]r]) = 1 : 3.7 (4) : 3.9 (4.3). (36) 

with BR{rjT]') = in the SU(3) flavor symmetry limit. The enhancement of the KK/tttt 
decay ratio relative to the naive SU(3) flavor estimate of Eq. (35) is a feature con- 
firmed by recent lattice results |2^. For the given choice of radial meson parameter 
{Rk = 2.82 GeV^^,Rrj = 2.48 GeV^^), as extracted from tensor meson decay, we 
obtain BR{KK) ^ BRirjr]). Choosing the SU(3) flavor limit of Rk = Rri leads to 



BR{KK) = 4.2 BR(riri), which is again consistent with the analysis of Ref. I^Qj. Hence, the 
proposed modelling of the Go decay mechanism contains the main features obtained from 
lattice calculations. 

As a measure of flavor symmetry breaking we introduce the parameter S defined by 
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S =< uu\V\g > / < ss\V\g >^ 1 leading to BR{riri') ^ 0. The flavor matrix elements for 
K K , rji] , rji]' decays listed in Table |I| correspond to the exact flavor symmetry limit {S = 1) 
and have to be generalised for pure glueball decays (aco = 1) 

M^„_,^^ = 2V2S 
M^^^^^ = V2 {als + S'Pls) (37) 
M^^_^^^, = 2 aps /3ps il-S'). 

If we adopt the extreme value of = 0.5 the strength of the KK and rji] decay 
modes will be suppressed by 50% and 34%, respectively, as compared to the unbroken flavor 
symmetry couplings. 

For scalar glueball masses at Mgq = 1300, 1500, 1700 MeV the reduced available phase 
space suppresses the rjr]' mode: 



Mgo [MeV] 


BR{7nv) 


: BR{KK) : 


BRiriri) : 


BR{r]i) 


1300 


1 : 


1.3 : 


2 : 


0.0 


1500 


1 : 


1.8 : 


2.6 : 


0.08 


1700 


1 : 


2 : 


2.9 : 


0.28. 



(38) 



Strong violation of flavor symmetry reduces the KK/titt ratio, while leading to a finite 
rjrj' contribution. However, data on charmonium decay involving hard gluons indicate no 
significant symmetry breaking [Q. We therefore choose to set S=l in the following discussion. 
Decay data on /o(1500) possess qualitative features arising both from the decay of a QQ 
(here dominantly nn) and a Go configuration, where a pure QQ or Go state interpretation of 
the /o(1500) is excluded. Given this simple first analysis of the decay pattern of the /o(1500) 
we are naturally led to a mixing scheme where both components are present. 

C. Decay analysis in the three-state mixing schemes 



In the following we give a decay analysis for the /o states defined by the three-state 
mixing schemes of Refs. [|] and p . Now, the couplings both of the QQ and Go components 
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to the two-meson decay channels are included. 

When the physical /o mesons are given as eigenstates of the interaction Hamiltonian Hj of 
Eq. (1), the partial decay widths of the process /o BC are defined by Eq. (5) with the 
ratio of coupling strengths k and the phase angle left as free parameters. The structure of 
Hi implies unbroken flavor symmetry, hence we choose S = 1. For the /o(1500) decays we 



give in Table |T| the squared decay amplitudes and the interference terms as defined in Eq. 
(5), where mass averaging as given by Eq. (6) is already included. The mass distribution of 
/o(1500) is parameterised with the values of Eq. (|26|) . 

The dependence of the partial decay widths T{BC) on the mixing coefficients {i = 
nn, ss, Gq) in the schemes of Ref. Q (= scheme A) and Ref. 0] (= scheme B) are studied by 
adjusting k and to the Crystal Barrel data for /o(1500) decays. First we have determined 
K and cos(f) from a fit to the observed decay modes /o(1500) KK and tttt. For both mixing 
schemes we obtain parameter values resulting in a /o — rjr] decay width in good agreement 
with the experimental number. This implies phenomenologically, that a Gq ^ GqGq contri- 
bution to the rjr] decay channel can be neglected. For this reason, we also take the rjr] decay 
mode into account when determining cos0 and k from a more restrictive fit. The partial de- 
cay widths have a rather weak dependence on the parameterisation of the mass distribution 
of the physical /o states. Therefore, the discussion on the detailed relationship between 
the /o(1370) and /o(400 — 1200) states is peripheral to our analysis. For our numerical 



results we use the resonance parameters M 



/o(1370) 



1360 MeK r 



/o(1370) 



350 MeV 



and M 



/o(mo) 



1697 MeV, F 



/o(1710) 



175 MeV fig. 



The mixing coefficients in scheme A [ffl are given as 



|/o(1370) > 
|/o(1500) > 
l/o(1710) > 



0.86 0.13 -0.5 
0.43 -0.61 0.61 
0.22 0.76 0.6 



\nn > 



\ss > 
\Go> 



(39) 



From a fit to the /o(1500) im- , K K , rjr] decay widths we deduce coscp = —0.67 and 
K = 0.095 GeV^^. Predictions for the partial decay widths of the physical /o states are 
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presented in Table For /o(1500) the predicted rjr]' decay width is considerably lower 
than the Crystal Barrel value, but the obtained ratio Tirjri') /T{riri) agrees with the GAMS 
value of 2.7 ± 0.8 for /o(1590) |ig. It is generally believed that /o(1500) and /o(1590) are 
manifestations of a single state, but then the large deviation in the experimental data for 
the rjrj' decay mode has to be clarified. Compared to the lowest order decay analysis, where 
only the coupling to the QQ component is retained, the pp mode is more suppressed and 
the An channel enhanced. We obtain BR{iin) : BR{pp) : BR{A'k) = 1 : 0.65 : 0.46 - 0.76 
for BR{'K* an) =0.1 — 1, still in deviation from the Crystal Barrel data. 
In the mixing scheme B the amplitudes for the physical states are given as 



|/o(1370) > 
|/o(1500) > 
l/o(1710) > 



0.84 0.28 -0.46 
0.40 -0.9 0.19 
0.36 0.34 0.87 



\nn > 



\ss > 
\Go> 



(40) 



From an analogous fit we obtain cos(f) = —0.953 and k, = 0.437 GeV^^. The predicted two- 
body partial decay widths are given in Table 0. Compared to the mixing scheme A, /o(1500) 
has a reduced Gq component, but the experimental data demand a significant Gq BC 
contribution in fitting the tttt, KK and 7]7] decays. This is balanced by giving the Go ^ BC 
transition a stronger weight than in the mixing scheme A, resulting in a higher value for k. 
Again, the predicted strength of the rji]' decay mode overestimates the experimental value, 
while similar results are obtained for the pp and Att decay modes. For mixing scheme B we 
obtain BR{mi) : BR{pp) : BR^An) = 1 : 0.46 : 0.55 - 0.97 for BR{n* an) = 0.1 - 1. 
In both schemes ain is a dominant decay mode with BR{nn) ^ BR{ain). With our fitting 
procedure, that is adjusting the relative strength of the Gq decay contribution to the vrvr, KK 
and rjf] decay channels of the /o(1500), similar results are obtained for the decay pattern 
of the /o(1500) in both mixing schemes. Particular predictions for the decay modes of the 
partner states /o(1370) and /o(1710) can be dramatically different in the two mixing schemes. 
In scheme A the total width for /o(1710) is in good agreement with the experimental result 
of r = 175 ± 9 MeV (F = 133 ± 14 MeV) [jlOl, but scheme B predicts a unobservable decay 
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width of r ^ 2000 MeV. Furthermore, for /o(1710) we get 

{0.93 for scheme A 
(41) 
0.60 for scheme B 

to be compared with the experimental value of 0.39 ± 0.14 ||10|| . 

The two mixing schemes predict rather similar results for the /o(1370) state. In both schemes 
/o(1370) has dominant pp and aivr decay modes. The ratio B R{K K) / {B R^rji]) is found to 
be 

{2.6 for scheme A 
(42) 
1.45 for scheme B 

to be compared with the Crystal Barrel value of B R{K K) / B R{'r]r]) ^2.5 
The main conclusions of this section remain unchanged if we use for /o(1500) a larger vrvr 
decay width as reported in Ref. [^], since a corresponding fit leads to similar values for k 
and cos(f) as for the adjustment to the Crystal Barrel data. 

Inclusion of the direct coupling of the Go component of the mixed /o(1500) state to the two- 
meson decay channels gives an improved description of its decay features. For both proposed 
mixing schemes predictions for the /o(1500) decays are very similar, since the difference in 
size of the Go amplitude is compensated by adjusting the strength of the Gq (QQ) (QQ) 
transition. 

Also, predictions for the /o(1370) state do not differ too much, since both schemes assign a 
dominant nn component to this state. The key difference rests on predictions for the partner 
state /o(1710). We have shown that a reduced Gq component in /o(1500), as proposed by 
Ref. 1^, is incompatible with the existence of a /q state at 1700 MeV as a narrow resonance. 
Conversely, if the existence of a narrow scalar component of the /j=o.2(1710) |T0[ state is 
confirmed, then a non- negligible glueball component is residing in the /o(1500) state. 
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D. Mixing amplitudes 



Given the shortcomings of the proposed mixing schemes of Refs. [^,^ we now directly 
extract the mixing amphtudes of the /q states in the given decay formahsm. We proceed 
in an analogous fashion, requiring a fit to the experimental data on /o(1500) decays, while 
demanding the existence of a narrow resonance /o(1710) as an additional requirement. 
First, we start with the unspecified mixing amplitudes of the /o(1500) state vector as defined 
by Eq. (2). The mixing amplitudes of /o(1500) are determined from a fit to the experimental 
data of the vrvr, KK and 7777 decay modes with the restriction of a weak pp decay channel. 
Exact fiavor symmetry [S = 1) for the quark-gluon coupling is assumed. Good agreement 
with the experimental data requires ann ■ dss < 0, i.e. a destructive interference between the 
nn and ss fiavor s, and we get 

|/o(1500) > = 0.314 \nn > - 0.581 \ss > + 0.751 |Go > (43) 

with the ratio of coupling strengths k = 0.1 GeV~^ and the phase factor cos(f) = —0.92. The 
fit to the /o(1500) decays alone does not fix the phase of the amplitude relative to the 
quarkonia coefficients. Instead we obtain the condition that -coscj) < 0. For the /o(1500) 
state vector we have taken > and coscj) < for reasons given later. The relative phases 
between the state amplitudes of Eq. (43) are the same as in mixing schemes A and B. For 
the full set of partial decay widths we obtain the results given in Table Compared to 
mixing scheme A the pp mode of /o(1500) is suppressed by a factor of two. This is achieved 
by enhancing the Go amplitude since the pp suppression for /o(1500) is driven by the gluonic 
component (see Sec. ^V]) . An alternative fit to obtain pp suppression requires a strong ss 
component and a larger value for k as resulting for example from the mixing scheme B. 
However, the existence of a narrow /o(1710) state rules out this possibility. 
For the A7t/2tt decay branching ratio we obtain B R{A7t) / B R(27i) = 0.57 — 1.04 using 
BR{tt* an) = 0.1 — 1, where recent results for tt* an favor the lower value. 
This is still in deviation from the Crystal Barrel value of BR{A'k) / BR{2'n) = 2.1 ± 0.6 [g. 
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The disagreement possibly hints at a sizable Go component in a. A strong coupling of cr to 
Go is supported by the experimental finding |T^ that aa is the largest hadronic branching 
ratio for the charmonium states xo and X2- In this case, the lowest order decay mechanism 
Go GqGq, neglected here, can enhance the aa and therefore the An decay mode. The pre- 
dicted ratio BR{riri') / BR{riri) = 2.2 is in agreement with the GAMS value of 2.7 ± 0.8 [|TU 



but in conflict with the Crystal Barrel result of 0.27 ±0.10 Breaking of flavor symmetry 
at the quark-gluon vertex, that is S" 7^ 1, will further enhance the tji]' decay width. 
With mass and eigenstate of /o(1500) flxed at 1500 MeV and by Eq. (43), respectively, we 
deduce for the bare masses rrii {i = nn, ss, Go) entering in Hj of Eq. (1) 

rrinn = 1500 MeV - 3.38 z, m,,- = 1500 MeV + 1.29 z, (44) 
mco = 1500 MeV + 0.18 z 

with the mixing strength 2; > 0. For the bare masses we obtain the level ordering rrinn < 
< iTiss, which is independent of z and corresponds to that proposed originally in Ref. [Q. 
When taking the alternative sign pattern of < and coscp > for the flt to the /o(1500) 
decays, we obtain the unphysical level ordering of the bare states with rrinn > rriss > t^Go-i 
we hence exclude this choice. Additionally, for the masses M<(>) of the mixing partners 
lying above (M>) and below (M<) /o(1500) we obtain: 

M> = 1500 MeV +2 z (45) 
M< = 1500 MeV - 3.9 z. 

For consistency, we call these states /o(1370) and /o(1710) and use the values for mass and 
total width given in the previous section for the parameterisation of the respective mass 
distributions. For the mixing partners of the /o(1500) we obtain the glueball-quarkonia 
content 

|/o(1710) > = 0.149 |nn> + 0.811 |ss > + 0.565 |Go > (46) 

and 
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|/o(1370) > = 0.938 |nn> + 0.070 |ss > - 0.341 |Go >, (47) 

where the amphtudes are independent of the particular value of z. Predictions for the cor- 
responding partial decay widths of /o(1710) and /o(1370) are also given in Table |V|. The 
total width and the ratio BR^titi)/ BR{KK) = 0.57 of the /o(1710) are in good agree- 
ment with the experimental results of of T = 175 ± 9 MeV (133 ± 14 MeV) |g and 
BR{'KTi)/BR{KK) = 0.39 ± 0.14 for /j=o,2(1710) For the /o(1370) the predicted ratio 
of BR{KK)/BR{riri) = 2.7 is to be compared to the experimental value ~ 2.5 p|. For 
the total width we get Ttot = 479 MeV, somewhat larger than the Crystal Barrel value of 
F = 351 ±41 MeV M- 



For a mixing energy 2; = 43 ± 31 MeV |T^, as deduced from lattice QCD in the quenched 
approximation, we obtain = 1508 ± 6 MeV, in agreement with the lattice result of 
tjigq = 1.65 ± 0.15 GeV 0. For the bare masses we get rrinn = 1355 ± 105 MeV and 
Mgs = 1556 ± 40 MeV. The values for the glueball-quarkonia mixing matrix element z fa- 



vored by previous approaches are z = 77 MeV 0, 2; = 64 ±13 MeV |T^ and z ^ 100 MeV, 
where later value is deduced from the full nonperturbative three-state mixing scheme of Ref. 
1^. Hence, semiphenomenological extractions of z favor a value near the upper limit set by 
the lattice result. 

Choosing an intermediate value of 2; ^ 80 MeV we obtain for the physical masses of the 
mixing partners given by Eq. (45), M> = 1660 MeV and M< = 1190 MeV, where M< 
is lower than the mass of the /o(1370). A similar result is obtained from the full mixing 
approach of Ref. [|], which corresponds qualitatively to our extracted scheme. 
Given our fit to the /o(1500) decays, we obtain a physical mass M< and a total width which 
deviate from the corresponding data on /o(1370). Therefore, within the framework of a 
three-state mixing scheme, the lower lying partner of /o(1500) could be identified with the 
very broad structure around 1100 MeV g], called /o(400 - 1200) Alternatively, the 

Crystal Barrel state /o(1370) could be the high mass tail of /o(400 — 1200). However, given 
the ansatz of Eq. (1), the most stable and testable consequences of our mixing scheme con- 
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cern predictions for all the two-body decay modes of the /o states as given in Table 0. Given 
the sensitivity of the decay pattern on the particular size of the Gq and QQ components in 
the /o states only, a full experimental determination provides a stringent test. 
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VI. SUMMARY AND CONCLUSIONS 



We have performed a detailed study of the two-body decay properties of the scalar- 
isoscalar /o(1370), /o(1500) and /o(1710) states resulting from the mixture of the lowest 
lying scalar glueball with the isoscalar quarkonia states of the 0''"'' nonet. In the decay 
analysis we have taken into account the coupling of the quarkonia and glueball components 
of the /o states to the quarkonia components of the two-meson final state, where the decay 
dynamics is guided by strong coupling QCD. Leading order corresponds to the transitions 
QQ {QQ) (QQ) and Go Go Gq, where the latter transition, which can contribute to 
rif], rji]' and aa final states, is omitted due to its ill-constrained nature. The coupling of the 
quarkonia component (QQ) of the /o states is described in the framework of the ^Pq pair 
creation model, which is successful in the phenomenology of OZI-allowed meson decay. In 
next-to-leading order we obtain a direct coupling of the Gq component of the /o states to the 
quarkonia component of the decay channel. The corresponding transition Go ~^ (QQ) (QQ) 
is modelled by resorting to a scalar digluonium wavefunction for Gq as given by cavity QCD. 
Predictions for the two-pseudoscalar decay channels, that is vrvr, KK and rji], are in rough 
agreement with recent lattice results in the limit of unbroken flavor SU(3). 
Our analysis on the two-body decay modes of the /o states focuses on the /o(1500), where 
detailed experimental data are available. In comparison to the observed decay properties 
of /o(1500) we elaborate the distinct decay patterns resulting from the unmixed QQ and 
Go components. Taking into account the coupling of the QQ components of the /o(1500) 
only, which corresponds to the leading order decay mechanism of Ref. f^, data on two- 
pseudoscalar decay channels are consistent with a dominant nn structure residing in this 
state. However, /o(1500) is too narrow for a simple nn structure, implying a sizable re- 
duction of its quarkonia components as resulting from glueball-quarkonia mixing. When 
going beyond the two-pseudoscalar decay channels large deviations from the observed de- 
cay pattern of the /o(1500) arise. The large Att branching ratio, fed by the aa and n*7r 
decay channels, and the suppression of pp is in clear conflict with the simple leading order 
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analysis of /o(1500), where coupling to its QQ component is taken, and points towards the 
relevance of a possible glueball component. Coupling of the Gq component of the /o(1500) 
to the quarkonia component of the two-body final state yields a significantly different decay 
pattern. In the two-pseudoscalar sector strong KK and rjrj decay modes are found. Further- 
more, 7r*7r and aa dominate over the vrvr decay channel, whereas pp is strongly suppressed. 
This latter feature of a strongly reduced pp coupling to Go is also obtained in a recent work 



31| , where the coupling of Go to the pionic decay channel is set up in an effective linear 
sigma model. Both observed features of the /o(1500), the enhancement of the 4-7?° decay 
channels and the reduction of the pp channel, are driven by its gluonic component. 
In a next step we have investigated the decay properties of the /o(1370), /o(1500) and 
/o(1710) states in the proposed three-state mixing schemes of Refs. [§|,|5[. Due to the inclu- 
sion of the Go {QQ) {QQ) mechanism in the decay analysis we were able to discriminate 
between the mixing schemes, in particular with respect to the orthogonal physical conse- 
quences for the /o(1710) state. A negligible Go component in /o(1500) as proposed by Ref. 
1^ is incompatible with the existence of a narrow /o state in the 1700 MeV mass region; 
hence our first analysis of the two-body decays favors the mixing scheme of Ref. [Q. 
Finally, we have extracted a three-state mixing scheme based on a detailed fit to the exper- 
imental /o(1500) decays, while demanding the existence of a narrow /o mixing partner in 
the 1700 MeV mass region. For the bare masses before mixing we obtain the level ordering 
rrinn < f^ss < TT^Goi Corresponding qualitatively to the scheme originally proposed in Ref. 0. 
Adjustment of the gluonic decay mechanism leads to a proper description of the observed 
two-pseudoscalar decay modes of /o(1500) with the exception of 7777'. Accordingly, we obtain 
a significant pp suppression coupled at the same time to the enhancement of the 4'7r''/27r 
ratio, in line with experimental requirements. A full reduction of the observed 47r/27r ratio 
requires a strong aa mode. Both deviations in the rjrj' and aa decay channels suggest a 
non-negligible leading order Go Go Go contribution, which was neglected here. A key 
signature for the extracted state composition of /o(1500) is the prediction for the oivr decay 
mode with BR{aiTT) ^ BR{TT7r). A pure QQ configuration or equivalently, application of the 
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leading order decay mechanism to the mixed /o(1500) state results in BR{ai7i) ^ 2BR{Tm). 
Hence an experimental determination of the aivr decay channel of /o(1500) is desirable to 
further quantify the relevance of its gluonic component. 

The resulting predictions for the decays of the /o(1370) and /o(1710) states are consistent 
with the limited experimental data P,pi]| on pseudoscalar decay modes. Clearly a search for 
the decay channels feeding the 47r decay products is most relevant in clarifying the nature of 
these partner states. Predictions for the total widths of /o(1370) and /o(1710) roughly agree 



with the observed values fgjl^]; the physical mass spectrum resulting from the extracted 
mixing scheme points towards a sizable influence of the broad /o structure at 1100 MeV [§, 
which should be clarified. 

The quantitative predictions of our analysis emphasise the role of the gluonic components 
in the /o states and its observable consequences in the two-body decay modes. Hence a full 
determination of the hadronic decay properties provides a sensitive test on the nature of the 
/o states, revealing the intrusion of the glueball ground state in the scalar-isoscalar meson 
spectrum. 
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APPENDIX A: SPATIAL WAVE FUNCTIONS 



For the spatial part of the wave functions used in the text we use harmonic oscillators 
^niiPhPj) which are defined in the limit of equal constituent masses as 



f f 1 

LT'\^{p^-myup^-P,) 



Pj. 



(Al) 



with the internal momentum and the size parameter R. The normalisation constant 
Nni depends on the radial and orbital quantum numbers n, 1 with 



nl 



1/2 



r{n + l + 3/2) 
The Laguerre polynomials are given by 



" t'o k\r{k + 1 + 3/2) ^ ■ 



(A2) 



(A3) 



For the QQ mesons we use: 



{nl) 



(00) for TT, rj, T]', K, K, p, uj 

(01) for oi, /o 
(10) for 7r*(1300) 



(A4) 
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APPENDIX B: QUARKONIA DECAY AMPLITUDES 



In the following we give the full expressions for the spin-spatial-color amplitudes 



M. 



{Ibc)SSC 
QQ-,BC 



defined in Eq. (12) as obtained for the quarkonia decay process QQ — > BC . 



For equal size parameters (= Rb) of the final state mesons B and C we obtain: 
1) QQ nn, rjr], rjrj', KK 



(Q)SSC 



QQ^BC 



I EbEc 



p5/2 p3 



B 



M 



fo 



7rV4 (i?2^ + 2i?|)V2 



exp 



^Bi^QQ + Rb) 



^5/2 




(Bl) 



2) QQ pp, ujuj 



M. 



{Ibc)SSC 



QQ-^BC 



M 



fo 



d5/2 d3 



exp 



QQ 

■2 



3 {^^m^ ' 




(B2) 



(B3) 



3) QQ 



aa 



M. 



{0)SSC 



I E„E„ 



p5/2 p5 
^QQ^B 



60K^/^ H ^ 



732 7d2 
^QQ^B 



RIq + 2i?l 



^ 



i^5/2_ 



RqqRURI 



Rl) 



(B4) 



;^9/2 



4) QQ 7r*(1300)7r 



QQ^7r*(1300)7r 



2'/' RTqRb 



p2 p2 



-^QQ + 2-R| 
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5) QQ oiTT 



+ 



p4 



^9/2 



(B5) 



(1)5SC 



exp 




\ K^'^ (B6) 
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APPENDIX C: TRANSFORMATION BRACKETS 

The relevant values for the transformation brackets < nl, NL, g\nili,n2l2, 9 > used in 
Eq. (19) are: 



1)^ = 



2)g^l 



< 00,10,0|01,01,0> = ^ (CI) 

v2 

< 01,01, 0|01, 01,0 > =0 (C2) 

< 10,00,0|01,01,0> = — ]= (C3) 

v2 



< 01,01,1|01,01,1 > = 1 (C4) 



< 00,02,2|01,01,2 > = ^ (C5) 

V2 

< 01,01, 2|01, 01,2 > =0 (C6) 

< 02,00,2|01,01,2 > = — ]= (C7) 

v2 
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APPENDIX D: DIGLUONIUM DECAY AMPLITUDES 



The spin-spatial-color matrix elements Mq^S^bc^ defined in Eq. (22), for the decay of 
the scalar digluonium state Gq (size parameter Rgq) into the final state mesons B C (size 
parameter Rb) are: 
1) Go ^ TTTT, 7777, W, KK (= 2 PS) 

?3/2 / p2 _ o d2 



2) Go ^ pp, (= 2 VM) 



3) Go — > (7(7 



4) Go ^ 7r*(1300)7r 



^Go-..*. - V 3^3 (i?| + 2i?^jV2 ^ ^^^^ 
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FIG. 1. Leading order decay mechanism of the Go — > GqGq transition. The dot indicates the 
three-gluon coupling. 

FIG. 2. Quark line diagram for the decay of the quarkonia component QQ into the final state 
mesons BC occurring as a leading order decay mechanism in strong coupling QCD. The dot indi- 
cates the ^Pq vertex with strength A. 

FIG. 3. Transition Go BC occurring as a next-to-leading order decay mechanism in strong 
coupling QCD. The dots indicate the quark-gluon coupling with strength 7. 

FIG. 4. Dependence of the branching ratios -Bi?(/o(1500) BC) on the mixing angle a, 
defined by the QQ configuration |/o(1500) >= cosa\nn > —sina\ss >, for the decay mechanism 
QQ BC of Fig. 2. For the meson size parameters we take the values given in the main text. 
The arrows indicate the predictions for the mixing schemes of Refs. in lowest order of the 
decay mechanism. 

FIG. 5. Dependence of the total width of the /o(1500) resonance on the mixing angle a defined 
by the QQ configuration |/o(1500) >= cosa\nn > —sina\ss >. The shadowed band corresponds 
to the experimental total width of F = 111 ± 12 MeV [|. 

FIG. 6. Dependence of the branching ratios BR{Cq BC) on the glueball mass for the decay 
mechanism Gq BC of Fig. 3. Exact fiavor SU (3) symmetry for the quark-gluon vertex is 
assumed. For the size parameters we use the values given in the main text. For mass averaging we 
take a glueball width of Fgq = 120 MeV. 



40 



TABLE I. Flavor matrix elements for the two-body decay QQ BC of the mixed fo{M) 
states as defined in Eq. (11). The quarkonia component of the /o states is defined as 
\QQ >= ann\nn > +ass\ss >. The quantities aps and Pps are related to the pseudoscalar 
mixing angle as defined by Eq. (13). For the vector mesons we use ideal mixing. 



TTTT, pp 



aa, LOLO 



7r*(1300)7r, aiTT \/6 Onn 

KK ann + a/2 ttss 

VV' 2 aps (3 PS {ann/\^ - ags) 
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TABLE II. Flavor matrix elements of Eq. (21) for the two-body decay Go — ^ BC of the mixed 

/o(M) states with glueball component ogqIGo >. The quark-gluon coupling is assumed to be flavor 
independent. 

B C M§„^sc/aGo 

TTTT, pp \/6 

(70", UU \/2 

7r*(1300)7r,ai7r VU 

KK 2 V2 

r/r/ \/2 

•qrf 
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TABLE III. Squared decay amplitudes and interference terms as defined by Eq. (5) for the 
transition /o(1500) — BC. Mass averaging as in Eq. (6) is already included. For the interference 
term we have Ibc = as restricted by the dynamical selection rule of the Go BC transition. 
All values are given in MeV. 



TTTT 52.2a2^ 224.4a|^ -103.1 an^^aco 

r)r) 20.9(0.63ann + 0.523asj)2 865.86^^, -131.1aGo(0.63ann + 0.523a5s) 

VV' 10.8{ann/V2~as~sf 

KK 10.7{ann + V^Gssf 828.0a^^ -86.8aGo{ann + V^a^s) 

PP 50.4a^^ IS.Oa^P 10.2annaGo 

uu 5.3a^^ 2.2a^g 2.0annaGo 

aiTT 94.2a2^ 

(7(7 264.6a^j, -63.1a„naGo 

7r*(1300)7r 7.2a2- 440.4a^ -55.2annaGo 
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TABLE IV. Predictions of mixing scheme A for the partial decay widths r{BC) of the 
physical /o states, B C are the final state mesons. The 47r decay channel refers to the con- 
tributions of the aa and 7r*7r intermediate states using BR{tt* — > crvr) = 0.1 — 1. All values 
are given in MeV. The predicted ratios Si?(/o(1370) K K) / BR{f o{l370) rjr]) = 2.6 and 
Si?(/o(1710) TT7r)/BR{fo{17W) KK) = 0.93 have to be compared to the experimental 
values of « 2.5 M and 0.39 ± 0.14 ITC 





^tot 


TTTT 


KK rjif] rjrj' 


PP 


LOUJ 


oivr 


* 

TT vr 


aa 


in 


Model 


136 


29.3 


/o(1500) 

3.9 5.2 19 


19 


2 


36.8 


9.9 


12.5 


13.5-22.4 


Exp. i 


111 ± 12 


32 ±8.4 


3.9 ±0.3 5.2 ±1.5 1.4 ±0.4 












68 ± 10.7 


Model 


410 


65.4 


/o(1370) 

29.1 10.8 3 


112.1 


22.8 


138.6 


12.7 


16 


17.3-28.7 


Exp. i 


351 ±41 




KK/r]ri « 2.5 














Model 


179 


23.4 


/o(1710) 

25.1 22.6 20.3 


21.3 


4.3 


28.2 


20.9 


12.7 


14.8-33.6 



Exp. lOl 133 ±14 vrvr/i^ri^: = 0.39 ±0.14 
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TABLE V. Predictions of mixing scheme B for the partial decay widths T{BC) 
of the physical /o states. The 47r decay channel refers to the contributions of the aa 
and 7r*7r intermediate states using BR{'n* — > (jtt) = 0.1 — 1. All values are given 
in MeV. The predicted ratios Si?(/o(1370) ^ K K)/BR{f o{1370) rjrj) = 1.45 and 

Si?(/o(1710) TT7r)/BR{fo{17W) KK) = 0.60 have to be compared to the experimental 
values of « 2.5 M and 0.39 ± 0.14 ITC 





^tot 


TTTT 


KK 


rjrj 


m' PP 


LOU) 


oivr 


* 

TT vr 


aa 




Model 


157 


34.7 


3.9 


5.1 


/o(1500) 

31.8 15.8 


1.6 


31.8 


16.2 


17.4 


19-33.6 


Exp. i 


111 ± 12 


32 ±8.4 


3.9 ±0.3 


5.2 ± 1.5 


1.4 ±0.4 










68 ± 10.7 


Model 


370 


38.7 


31.4 


21.6 


/o(1370) 

1.3 110.8 


22.6 


132.2 


8.2 


2.8 


3.6-11 


Exp. i 


351 ±41 




KK/iji] 


K, 2.5 














Model 


2000 


253 


381 


442 


/o(1710) 

0.4 63 


12 


76 


536 


263 


316-799 



Exp. lOl 133 ±14 vrvr/i^ri^: = 0.39 ±0.14 
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TABLE VI. Predictions for the partial decay widths T{BC) of the physical /o states 
for the extracted mixing scheme. The 47r decay channel refers to the contributions of 
the aa and 7r*7r intermediate states using BR{tt* — > (Ttt) = 0.1 — 1. All values are 
given in MeV. The predicted ratios BR{ fo{l370) KK)/BR{fo{l370) r]r]) = 2.7 and 
Si?(/o(1710) TT7r)/BR{fo{17W) KK) = 0.57 have to be compared to the experimental 
values of « 2.5 M and 0.39 ± 0.14 ITC 





^tot 


TTTT 


KK 


rjif] 


m' PP 


LOU) 


oivr 


* 

TT vr 






Model 


101 


23 


2.8 


6.7 


/o(1500) 

14.6 9.7 


1 


19.6 


11.8 


12.1 


13.3-23.9 


Exp. i 


111 ± 12 


32 ±8.4 


3.9 ±0.3 


5.2 ± 1.5 


1.4 ±0.4 










68 ± 10.7 


Model 


479 


78.2 


27.4 


10.1 


/o(1370) 

4.6 133.2 


27.1 


164.9 


14.6 


19 


20.5-33.6 


Exp. i 


351 ±41 




KK/iji] 


K, 2.5 














Model 


146 


15 


26.3 


24.5 


/o(1710) 

27.7 9.7 


1.9 


13 


18 


10.1 


11.9-28.1 



Exp. lOl 133 ±14 vrvr/i^ri^: = 0.39 ±0.14 
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